Measuring Lagrangian accelerations using an instrumented particle 
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Accessing and characterizing a flow impose a number of constraints on the employed measurement 
techniques; in particular optical methods require transparent fluids and windows in the vessel. 
Whereas one can adapt apparatus, fluid and methods in the lab to these constraints, this is hardly 
possible for industrial mixers. We present in this article a novel measurement technique which is 
suitable for opaque or granular flows: an instrumented particle, which continuously transmits the 
force/acceleration acting on it as it is advected in a flow. Its density is adjustable for a wide range of 
fluids and because of its small size and its wireless data transmission, the system can be used both in 
industrial and scientific mixers allowing a better understanding of the flow within. We demonstrate 
the capabilities and precision of the particle by comparing its transmitted acceleration to alternative 
measurements, in particular in the case of a turbulent von Karman flow. Our technique shows to 
be an efficient and fast tool to characterize flows. 



I. INTRODUCTION 

Experimental fluid dynamics research in the lab 
consists of an interplay of suitable flow generation 
devices, working fluids, measurement techniques 
and analysis, with goals ranging from fundamental 
research in statistical / non-linear physics to the 
optimization of mixers in industrial R&D depart- 
ments. In this endeavor, very significant progress 
has been achieved during the last decade with the 
advent of space and time resolved optical tech- 
niques based on high speed imaging [1]. However, 
direct imaging is not always possible especially in 
industry: opaque vessels, non-transparent fluids, 
environmental constraints among other things may 
be limiting factors. Even if the fluid is transparent, 
the injection of tracer particles might be still not 
allowed or unsuitable due to bio-medical or food 
regulations, or due to the chemical properties of 
the fluid. While techniques using other kinds of 
probing waves {e.g. acoustics [2]) have been de- 
veloped, a direct resolution of the Eulerian flow 
pattern is not always possible. In this context, 
Lagrangian techniques provide an interesting al- 
ternative particularly for problems related to mix- 
ing [3, 4]. 

Lagrangian tracers with a temperature sensi- 
tive dependance have been used in the study 
of Rayleigh-Benard convection [5], a problem for 
which our group has developed instrumented par- 
ticles [6-9]. The approach was to instrument a 
neutrally buoyant particle in such a way that it 
measures the temperature fluctuations during its 
motion as it is entrained by the flow, while trans- 
mitting the data via radio frequency to a lab opera- 
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tor in real time. Meaningful information regarding 
the statistics of thermion plumes have been ob- 
tained, with excellent agreement with other tech- 
niques [5] and direct numerical simulations [10]. 
In the work reported here, we built upon this ap- 
proach to instrument the particle such that one 
gets flow parameters directly from the measure- 
ments (in [6], one had to simultaneously film the 
particle motion). We equip the particle with a 3- 
axis accelerometer, whose measurements are sam- 
pled at a rate equal to 316 Hz and transmitted to 
the lab operator. This particle is intended for tur- 
bulent flows. Thanks to its radio transmission it is 
suitable for opaque fluids or apparatuses without 
access for optical measurement techniques. Its con- 
tinuous operation is also advantageous over Par- 
ticle Tracking Techniques which have to operate 
in chunks as the memory of the tracking cameras 
is necessarily limited. Moreover and in contrast 
to tracer particles this instrumented particle can 
be easily re-extracted from the apparatus after the 
experiment. However, as the particle is advected 
in a flow it rotates and consequently continuously 
changes its orientation with respect to the labo- 
ratory frame. Thereby the signals of the 3D ac- 
celerometer are altered in a non-trivial way, and 
detailed characterization and methods to extract 
meaningful information from the acceleration sig- 
nals are needed. We present here the preliminary 
results of this characterization. 

This article is organized as follows: first, we 
present the instrumented particle and additional 
techniques needed for its characterization (sec- 
tion II). In section III, we present an analysis of 
the results obtained in two different configurations: 
First, a simple pendulum with the particle at- 
tached at the end of a stiff arm, then the particle 
advected in a fully turbulent flow. In order to ver- 
ify that the transmitted acceleration is well related 
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to its motion, we compare the results to simultane- 
ous alternative measurements. Finally, we discuss 
limitations and perspectives of this new measure- 
ment technique (section IV). 

II. "SMART PARTICLES" 

The apparatus described in the following is de- 
signed and built by smartlNST S.A.S., a young 
startup situated on the ENS de Lyon campus. The 
device consists of a spherical particle (the so-called 
smartPART) which embarks an autonomous cir- 
cuit with SD-acceleration sensor, a coin cell and a 
wireless transmission system; and a data acquisi- 
tion center (the so-called smartCENTER), which 
acquires, decodes, processes and stores the signal 
of the smartPART (see Fig. 1). The ensemble - 
smartPART and smartCENTER - measures, dis- 
plays and stores the three dimensional acceleration 
vectors acting on the particle as it is advected in 
the flow. The accelerations are observed in a mov- 
ing and rotating coordinate system and consist of 
four contributions: gravity, translation, noise and 
possibly a weak contribution of the rotation around 
the center of the particle itself. 



A. Design Technical Details 

a. Sensor: The central component of the 
particle is the ADXL 330 (Analog Device) - a 
three axis accelerometer. This component be- 
longs to the category of micro-electro-mechanical 
systems (MEMS). Each of the three axes returns 
a voltage proportional to the force acting on a 
small, movably mounted mass-load suspended 
by micro- fabricated springs. The three axes 
of the ADXL 330 are decoupled and form an 
orthogonal coordinate system attached to the 
chip package. From this construction arises a 
permanent measurement of the gravitational 
force/acceleration g = 9.8m/s^ • = g • Cg. Each 
axis has a guaranteed minimum full-scale range 
of ±3 ^; however, we observe a typical range of 
±3.6^ = 35m/s^ per axis. The sensor has to be 
calibrated to compute the physical accelerations 
from the voltages of the accelerometer. 

b. smartPART: The signals from the ADXL 
330 are first-order low-pass filtered at fc = 160 Hz 
and then digitized at 12 bits and 316 Hz sampling 
rate. A multiplexer prior the signal digitization 
induces a small time shift between the compo- 
nents of 0.64 ms. The output is then reshaped 
into small packets and send via radio frequency. 
The ensemble is powered by a coin cell. A voltage 



regulator ensures a stable supply voltage and thus 
a constant quality of the measurement. A Hall 
switch allows one to power-down most compo- 
nents; thereby, the battery is only used during 
experiments. Depending on the power needed to 
transmit the acceleration signals, a particle oper- 
ates continuously for 6 to 36 hours. The ADXL 
330 is soldered to the printed circuit board such 
that it is situated close to the geometrical center 
of the particle. The particle itself is spherical 
with a diameter of 25 mm. The capsule walls are 
made of Polyether-ether-ketone (PEEK) which is 
known for its excellent mechanical and chemical 
robustness. It is leak-proof and its density can be 
matched to fluids by adding extra weight (namely 
Tungsten paste) to the particle's interior; within 
the density range of 0.8 — 1.4 g/cm^ a relative 
density match of better than 10~^ is achievable. 
The particle is thus suited for most experiments 
in water and water-based solutions. It should 
be noted that the mass distribution inside the 
particle is neither homogeneous nor isotropic: in 
particular its center of mass does not coincide with 
the geometrical center, making it out-of-balance. 
In practice this results into a pendulum-like 
motion of the particle in the flow. Nevertheless, 
the imbalance can be adjusted to some extent by 
adding patches of Tungsten paste to its interior, 
and the particles we use are carefully prepared 
such that they are neutrally buoyant, avoid any 
pendulum-like behavior and rotate easily in the 
flow. 



c. smartCENTER: The signals from the 
smartPART are received by an antenna connected 
to the smartCENTER, which contains radio re- 
ception, processing and display units. It demod- 
ulates and decodes in real-time the received raw 
signal into a time-series of raw voltages of the 
ADXL 330. The physical acceleration sensed by 
the smartPART asp can then be computed: 

/ ai \ / (Ai - Oi)/Si \ 
asp = U2 = {A2 - 02)/S2 , (1) 
\a^l V(^3-03)/^3 / 

where A^, Oi and Si are the measured raw signal, 
the offset and the sensitivity of each axis, respec- 
tively. Offset and sensitivity have to be calibrated 
beforehand; the procedure is described in the fol- 
lowing section. The resulting time-series are saved 
for further processing. 
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Figure 1: a) and b) : the instrumented particle (smartPART) and its data control & acquisition unit 
(smartCENTER) . The coin cell is 20mm in diameter. c) The diagram sketches the measurement of the 
acceleration, its transmission to and the processing by the smartCenter. 



Calibration and robustness 



The offset and sensitivity of the ADXL 330 have 
to be calibrated to convert the measured voltages 
into a physical acceleration. The axes of the ac- 
celerometer form an orthogonal coordinate system 
according to Eq. (1). At rest one observes only 
gravity projected onto the sensor at an arbitrary 
orientation. The observed raw values define con- 
sequently a translated ellipsoid (for simplicity we 
set 1^1 = 1): 



We find that the particle at rest has an average 
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Eq. (2) can be arranged to: 



(2) 



(3) 



with six parameters containing offset and sen- 
sitivity. A sufficient number of measurements 
with different orientations define a set of equations 
which is solved using a linear least squares tech- 
nique. Offset and sensitivity are then 
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and Sj = 
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noise of (Jx = ay = 0.006 g and az = 0.008 giving 
1^1 = ~ 0-012 g^ where g is the magnitude 

of gravity. An analysis using the residuals showed 
a slightly higher resolution of (Jx = (Jy = 0.005 g 
and az = 0.003^, and \cr\ = 0.008^. These values 
are thus the absolute errors of our measurement. 

The ADXL 300 has among other things been 
chosen for its weak temperature dependance: its 
offset typically varies by 10-^ ^/°C, and its sensi- 
tivity by 0.015 %/°C. Digitizing and transmission 
unit were verified to be temperature independent. 
Consequently, the total temperature dependence 
of the smartPART is given by its accelerometer. 
For high precision measurements, it is advised to 
calibrate the particle at experiment temperature 
shortly before the actual experiment. 

We noticed a small drift of the order of 0.005 g/h 
for the 2:— axis. No drift was observed for the x— 
and y— ax.es. Since a voltage regulator ensures a 
stable supply voltage of the circuit, this drift stems 
most likely from the internal construction of the ac- 
celerometer. Owing to the continuous data trans- 
mission of the instrumented particle, one flow con- 
figuration can be characterized in approximately 
30 minutes. Hence, the little drift of the axis 
can be neglected. 

Considering the mechanical robustness, the 
smartPART survived several days in a von Kar- 
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man mixer and neither contacts with the wah nor 
with the sharp edged blades of the fast rotating 
propehers damaged its function or sheh. 



III. ACCELERATION SIGNALS 

As mentioned previously, the smartPART trans- 
mits in real-time the accelerations acting on the 
particle as it is advected in the flow. The noise-to- 
signal ratio being small, we neglect the noise from 
here on. The contributions consist therefore of: 
gravity, translation, and rotation of the particle. 
We now test the accuracy of the particle signals in 
two different experimental configurations by com- 
parison with alternative measurements. 



A. 2D pendulum 

A pendulum is a simple and well-known case, 
ideal to measure the resolution of the particle. 
A stiff pendulum with 60 cm long stiff arm is 
equipped with a position sensor returning the 
deflection angle ip of the arm. The particle is 
fixed at known length, with a known arbitrary 
orientation to the arm. The fact that a rotation of 
the particle around its center is restricted, implies 
that the contribution from the rotation of the 
particle around is center is known. Measuring asp 
at rest ((/p = 0°) and at several arbitrary positions 
one can determine the axis of rotation of the arm. 
Once this vector is known the measured accel- 
eration signal is rotated/re-expressed such that 
ay points with the arm, ax with the movement, 
and az with the axis of rotation. Note that by 
definition the latter does not change when the 
pendulum moves. The signal seen by the particle 
is a two-dimensional problem and fully described 
as a function of the deflection angle if. 
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Figure 2: Comparison of the particle's signal asp(t) 

(•) to the theoretical curves based on the position 
sensor ( — ). ay points with the arm and ax measures 
the force in direction of the movement. No force is 
exerted along the ^— axis (the green lines represent 
the uncertainty of the calibration). Note, that 
acceleration is measured in ^ = 9.8m/s^ . 

Fig. 2 shows asp(t) for several periods of the 
pendulum, measured by the smartPART and by 
the position sensor. The agreement between the 
two signals is very good, and in particular better 
than the uncertainty of the calibration. Hence, the 
Lagrangian acceleration of the smartPART corre- 
sponds well to its actual motion in this simple case. 
We now move on to more complicated motions. 

B. Fully developed (3D) turbulence 



(ix{^) =9 sin (/:) + 
ay{Lp) = g cos if -\- lif"^ . 

In the limit of the small angle approximation, 
this simplifies to the well-known oscillations of fre- 
quency u: 

(^x{^) ^ — /cj^ sincjt, 

ay [(f) ^ g -\ — ^ V 1 + 2ujt) . 

The simultaneous measurement of the angle, (p{t)^ 
and the particle's signal, asp(t), enables us to com- 
pare the two signals without any other approxima- 
tion or fit than Eq. (5). 



The instrumented particle is intended for the 
characterization of complex/turbulent flows. Such 
flows exhibit strong, intermittent variations in the 
acceleration. To verify the suitability of the smart- 
PART for these conditions, we now investigate its 
motion in a fully turbulent mixer while tracking it 
with an independent optical technique. 

Namely, we use a von Karman water flow: a 
swirling flow is created in a square tank by two 
opposing counter-rotating impellers of radius R = 
9.5 cm fitted with straight blades 1cm in height. 
The flow domain in between the impellers has char- 
acteristic length H = 20 cm = 2R (see Fig. 3) and 
the vessel is built with transparent flat side walls, 
allowing direct optical measurements over almost 
the whole flow domain. Blades on the impellers 
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Figure 3: Experimental setup: a) picture of the apparatus; b) sketch of the arrangement; c) a textured 

instrumented particle at different orientations. 




Figure 4: A sample trajectory of the instrumented particle seen by the camera ( — ) or smartPART (•), it is 
/prop = 3 Hz. The absolute orientation enables us to re-express the camera measurement of the particle (lab 
frame) in the moving frame of the particle and vice-versa. In the former gravity is subtracted and in the particle 

frame gravity is represented by the red line. 
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work similar to a centrifugal pump and add a 
poloidal circulation at each impeller. For counter- 
rotating impellers, this type of flow is known to 
exhibit fully developed turbulence [11]. Within a 
small region in the center the mean flow is little 
and the local characteristics approximate homoge- 
neous turbulence. However, at a large scale it is 
known to have a large scale anisotropy [12, 13]. 
At a propeller frequency of 3 Hz we estimate a 
Reynolds number based on the Taylor micro-scale 
of Rx = 500 ± 50. 

We optically track the translation and absolute 
orientation of the smartPART while simultane- 
ously acquiring the transmitted acceleration time- 
series. These optical measurements are then used 
as a reference to compare with the instrumented 
particle's signals. The six-dimensional tracking 
technique (or 6D tracking, 3 components for the 
translation and 3 components for the rotation of 
the particle around its center) is explained in de- 
tail in [14, 15] and briefly sketched here (Fig. 3). 
In order to determine the absolute orientation the 
particle is textured by hand using black-ink per- 
manent marker (see Fig. 3c). Acceleration sen- 
sor and texture are then calibrated/retrieved in- 
dependently; nevertheless, the accelerometer is at 
a fixed but unknown orientation with respect to 
this texture, i.e. sensor and texture are related 
by a constant rotation matrix. We determine this 
matrix by acquiring the acceleration signals of the 
particle at arbitrary orientations while addition- 
ally determining its orientation and the location 
of gravity on the texture. The particle is then in- 
serted into the apparatus, which is illuminated by 
high power LEDs. Its motion is tracked by two 
high-speed video cameras (Phantom VI 2, Vision 
Research) which record synchronously two views 
at approximately 90 degrees. The observation vol- 
ume isl5xl5xl5 [cm^] in size and resolved at 
a resolution of 4.2 pixel/mm. In our configuration, 
a camera can store on the order of 14 000 frames 
in on-board memory, thus limiting the duration of 
continuous tracks. Therefore, a computer issues 
the recording of 8 bit gray-scale movies at a suffi- 
ciently high frame rate while controlling the smart- 
CENTER such that acceleration signal and images 
are synchronized. After extracting the time-series 
of the particle's position and orientation, one can 
then compare the accelerometer's signal to the mo- 
tion of the particle. 

It should be stressed that the two measure- 
ment techniques observe the motion of the instru- 
mented particle in two completely different refer- 
ence frames. On the one hand, the 6D tracking 
uses a fixed, non-rotating coordinate system, and 
is referred to as the lab frame. On the other hand, 
as the particle is advected and turned in the flow. 



it and consequently the embarked accelerometer 
constantly rotate their coordinate system with re- 
spect to the lab frame; the acceleration signal is 
thus measured in a frame which is continuously 
rotating and not fixed. This frame is referred to as 
the particle frame. The acceleration sensor mea- 
sures the forces acting on it as it moves in the flow. 
Knowing the absolute orientation of the particle at 
each instant we can express the signal of the smart- 
PART in the lab frame by rotating it such that it 
corresponds to a non-rotating particle. Starting 
from the time-series of position and orientation, 
it is also possible to compute the linear, centrifu- 
gal and gravitational acceleration/force acting on a 
point inside the particle and then project these into 
the rotating particle frame. The different compo- 
nents are then expressed in the frame of the sensor. 

Fig. 4 shows a sample trajectory in both coor- 
dinate systems. The agreement between the two 
techniques is remarkable. Furthermore, one ob- 
serves that the projection of gravity is continuously 
changing: the particle is rotating in a non-trivial 
way. Deviations between the two techniques stem 
from several experimental errors. First, the posi- 
tion measurement: bubbles, reflections and other 
impurities alter the measured position of the par- 
ticle. The acceleration is the second derivate and 
thus highly sensitive to such events. Second, the 
orientation measurement: the absolute orientation 
is needed to change between the reference frames. 
The uncertainty in the absolute orientation is typ- 
ically 3°; that results in a wrong projection of 
gravity of ±0.5m/s^. It further biases the rota- 
tional forces, as they are derivatives of the ori- 
entation time-series. Finally, the matrix relating 
sensor and texture: this matrix is constant and 
thus a systematic contribution. The uncertainty 
is less than 2° - i.e. the error in projecting grav- 
ity is < 0.3 m/s^. The observed agreement in the 
lab frame, Aa = asp — aeo, between the two 
techniques is as follows: all three components of 
Aa have the same PDF. Surprisingly, the (ab- 
solute) uncertainty almost doubles by increasing 
/prop from 2 Hz to 3 Hz. Nevertheless, for 80% 
of the data the agreement is better than 0.8 m/s^ 
and 1.6 m/s^, respectively. For comparison, the 
absolute value laeol has a mean of 2.9 m/s^ and 
6.6 m/s^ and a standard deviation of 1.8 m/s^ and 
4.1 m/s^, respectively. The signal of the particle 
is thus corresponding to the flow, however, its in- 
terpretation is not simple. In particular, and after 
comparing many different trajectories, it becomes 
clear that no easy transformation is available to 
get rid of the rotation of the particle. 

By construction the center of the accelerometer 
is placed at r = 3 mm • e^. A rotation of the par- 
ticle around its geometric center will thus add a 
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Figure 5: Ratio of the rotational forces to the total Figure 6: Preferred position of the instrumented 

force acting on the particle. The 80% percentile is particle: independent of the propeller speed it is 

found at a ratio of 0.14 and 0.16, respectively. mostly situated in a torus shape around the propeller. 



centrifugal contribution to the measured accelera- 
tion. This rises the question which term - trans- 
lation or rotation of the particle - dominates the 
acceleration signal. To address this question we 
take advantage of the 6D-tracking, which enables 
us to compute the different forces acting on a point 
at r = 3 mm • inside the sphere. We can thus 
compare the contribution of the translation and 
that of the rotation of the particle. Fig. 5 shows 
the ratio of the rotational {i.e. centrifugal) accel- 
eration, arot = u^xujxr-\-^xr to the total 
acceleration, atrans + o^rot, (without gravity). Di- 
mensional arguments tell that atrans /prop and 
<^rot /prop- Consistently, the PDF of the ratio 
I o^rot I / 1 citrans + cirot | differs ouly little for the two 
propeller frequencies. Moreover, it is peaked at 
5% and the 80% percentile is at a ratio of 14% and 
16%, respectively. Hence, it is legitimate to neglect 
the rotational forces if no 6D tracking is available. 



IV. DISCUSSION 

In the latter part of this article we studied the 
behavior of a large neutrally buoyant sphere in a 
turbulent flow. Comparing with solid spheres of 
the same size in the same mixer, we find that the 
particle in general behaves almost identically [16]. 
In particular (and despite the fact that the in- 
strumented particle is neutrally buoyant) we ob- 
serve that it generally stays in a region close to 
the impellers. Fig. 6 shows the PDF of position 
for the smartPART. Independent of the impeller 
speed it is mostly situated in a torus shape around 
the propeller, exhibiting a preferential sampling of 
the flow for these large neutrally buoyant spheres. 

Moreover, since we investigate large particles 
with a size Impart comparable to the integral length 
scale. Lint, moving through the whole mixer, the 



Kolmogorov assumptions to characterize turbu- 
lence are no longer valid. For these reasons, the 
smartPART can be insufficient to access all details 
of a turbulent flow: some parts of the flow are little 
explored, and some scales of the turbulence might 
be filtered due to the size of the instrumented par- 
ticle. However, one should bear in mind that these 
features of the flow are often accessed by means of 
optical methods whereas the instrumented particle 
operates also in environments and fluids which are 
unsuitable for optical measurement techniques. 

Some other experimental constraints should be 
additionally stressed here. As previously said, the 
mass distribution inside the particle is neither ho- 
mogeneous nor isotropic. It is therefore possible 
that the particle is out-of-balance, i.e. that the 
center of mass does not coincide with the geomet- 
rical center. Such a particle has a strong preferred 
orientation and wobbles similar to a kicked physi- 
cal pendulum. The imbalance can be adjusted to 
some extent by adding weight to its interior, but 
the particle must be prepared very carefully and 
one must make sure that the particle used is well- 
balanced and rotates easily in the flow. 

Also, the receiver/demodulation unit of the 
smartCENTER works best within a range of 
radio power, i.e. particles which are emitting 
either too strong or too weak are undesirable 
and one has to adjust the radio emission of the 
smartPART. A stronger radio emission power 
can be required, e.g. if the apparatus builds 
a Faraday cage {i.e. an electrically-connected 
metal structure surrounds the flow), or if the 
signal has to pass a longer distance in more 
water or in a bigger apparatus. Solutions with 
a high conductivity are also likely to damp the 
radio signal. Naturally, a stronger radio emission 
shortens the life time of the battery. Neverthe- 
less, particles with stronger radio emission still 
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last 6 to 12 hours, which is sufficient in most cases. 

To conclude, we presented the working principle 
of an instrumented particle giving a measure of 
the three components of the Lagrangian accelera- 
tion. We were able to show that the Lagrangian 
acceleration of the smartPART corresponds well 
to its actual translation and is not biased by a 
possible rotation of the particle around its center. 
Work on extracting detailed information on the 
flow from the acceleration time-series is ongoing. 
This instrumented particles can shed some light 
into mixers which were not or hardly accessible 
up to now. Due to its continuous transmission 
one flow conflguration can be characterized within 
~ 30min. Apart from its appeal for chemical 
and pharmaceutical industry, it might be an 
interesting tool to quantify flows in labs. 
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